Quasi-deterministic generation of maximally entangled states of 
two mesoscopic atomic ensembles by adiabatic quantum feedback 
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We introduce an efficient, quasi-deterministic scheme to generate maximally entangled states of 
two atomic ensembles. The scheme is based on quantum non-demolition measurements of total 
atomic populations and on adiabatic quantum feedback conditioned by the measurements outputs. 
The high efficiency of the scheme is tested and confirmed numerically for ideal photo-detection as 
well as in the presence of losses. 
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I. INTRODUCTION 

Quantum entanglement is one of the most fundamen- 
tal aspects of quantum mechanics, as well as an essen- 
tial resource in quantum communication and informa- 
tion processing. Although very difficult to realize, en- 
tangled states of material particles have been thoroughly 
studied in recent years both theoretically and experimen- 
tally, and some schemes for their generation have been 
designed and partially realized. Some studies concen- 
trated on how to produce entanglement between groups 
of two or few atoms, exploiting for example the collec- 
tive vibrational motion of trapped ions Jl|, |2( , the single- 
photon interference at photodetectors @], or the condi- 
tional dynamics of two atoms within a single-mode cavity 
field Q. More recently, there has been a growing inter- 
est on how to create multipartite entanglement between 
atoms belonging to a single atomic ensemble considered 
as a multi-party quantum system, by exploiting the inter- 
action with a light field [f| , and the subsequent detection 
process 0, Q- Finally, and more ambitiously, various 
schemes have been proposed for the entanglement of dif- 
ferent (two or more) macroscopic or mesoscopic atomic 
ensembles |1 S EH GH El 

In the cases of several (at least and typically two) 
macroscopic atomic ensembles, where collective atomic 
operators can be described by some continuous-variable 
approximation, it is only possible to design schemes for 
the realization of weak entangled states. Some of these 
schemes exploit quantum non-demolition (QND) mea- 
surements on auxiliary electromagnetic fields (usually 
assumed in some Gaussian state) interacting with the 
atoms to prepare entangled states of atomic systems 
[EEHH} However, the probabilistic nature of the quan- 
tum measurement events makes the generation of atomic 
entangled state conditioned by the measurement out- 
comes, usually yielding a low probability of success. This 
is particularly true for the preparation of maximally en- 
tangled states. Such a shortcoming should be in principle 
overcome by exploiting the knowledge of the state vec- 
tor of the atomic system conditioned on the outcome of 
a measurement, and then by introducing a proper feed- 



back scheme to efficiently drive the system toward a max- 
imally entangled state. Actually, Stockton et ai. showed 
in Ref. Q how this strategy can be properly used to de- 
terministically prepare highly entangled Dicke states of 
a single atomic ensemble. 

In the present work we introduce a reliable feedback 
scheme to generate maximal entanglement of two meso- 
scopic atomic ensembles. In this scheme, the discrete 
quantum nature of the atomic systems is fully taken into 
account without resorting to any continuous variable or 
Gaussian approximation. Our proposal is based on the 
model introduced by Di Lisi and M0lmer E2j where 
two collections of atoms, probed by a sequence of single- 
photon scattering processes, are conditionally entangled 
by QND measurements of the total atomic population 
difference between the two atomic samples. This model 
has been recently shown to be robust against spontaneous 
scattering Ell- 
in the present work the results of the QND measure- 
ments obtained by photo-detections are exploited to drive 
the system into the maximally entangled state by a suit- 
able feedback mechanism. The feedback scheme that we 
introduce is a proper modification to the fully discrete 
case of the continuous feedback strategy originally de- 
signed by Thomsen, Mancini and Wiseman (TMW) EH) 
to generate high spin squeezing of a single atomic ensem- 
ble, whose experimental realization was recently obtained 
by Geremia et aZ|l5l|. The same scheme was generalized 
to the case of two atomic ensembles to produce two-mode 
spin squeezing |T(| . 

Our procedure is monitored by quantitative wave- 
function simulations which show how the sequence of 
photo-detection events, followed by the feedback signal, 
gradually modifies the state of the samples and post- 
selects the maximally entangled states. We show that 
the feedback scheme enormously increases the rate of suc- 
cess in producing maximally entangled states of the two 
atomic ensembles compared with the scheme in which 
feedback is absent. We also show that the efficiency is 
further improved by adiabatically switching off the feed- 
back signal; in this way one obtains a quasi-deterministic 
generation of the maximally entangled state. Finally, we 
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FIG. 1: (Color online). Level structure of the atoms. \g) 
and |/) are metastable states, \g) is coupled off resonantly by 
the electromagnetic field, whose annihilation operator is ai, 
to the excited state \e). Here 7 is the spontaneous transition 
rate and A is the detuning between the coupling field and the 
atomic transition frequency. 



study the problem for the more realistic case of imperfect 
detectors, and we show how the feedback scheme guar- 
antees a very high probability of success in this case as 
well, making the mechanism quite reliable against losses. 




FIG. 2: (Color online). Schematic experimental configura- 
tion. Atoms occupying the internal state \g) in the two sam- 
ples interact with the light field, incident from the left. The 
presence of the cavities guarantees a high spectral and direc- 
tional resolution of the entering photon. The phase shift of 
the light field due to the interaction with the atoms is reg- 
istered by the different photocurrents in the two detectors. 
These signals can be combined with that produced by a sig- 
nal generator and hence sent back to the two ensembles. 



II. THE MODEL 

The two atomic ensembles, denoted by "1" and "2" 
respectively, are identical, and each one is constituted by 
N identical atoms in a static magnetic field, whose level 
structure consists of two metastable lower states, \g) and 
I/), that correspond to Zeeman sublevels of the electronic 
ground state of alkali atoms, and one excited state |e), 
cf. FigHI 

We then consider an optical beam passing through the 
atomic clouds which is coupled (out of resonance, with 
detuning A) only to the transition \g) — > |e). We can 
introduce the atomic spin operators for an atom a in 
ensemble i as 
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where a = 1, ...,N is the atomic index and i = 1,2 is 
the ensemble index. The dynamics of the ensembles can 
be described by collective spin operators whose x-,y- and 
z-components, for each ensemble i, read 
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These sums define the x-,y- and z-components of the col- 
lective angular momentum jW. In particular, the eigen- 

*(i) 

values rrii of Jz are proportional to the population dif- 
ference in the two stable states. 

The two atomic ensembles are initially prepared, by 
optical pumping, so that each ensemble is fully polarized 

al ong the x-axis with collective spin equal to = N/2 
0[iilD3[iJ This means that the atoms arc distributed 
between states \g) and |/) according to a binomial distri- 
bution with probability 1 /2 for each state. The composed 
system made by the two atomic clouds is described by the 
total spin operators 



J± = ) ± jf 5 , with k 



x,y,z . 



(3) 



The maximally entangled state of the composed system 
can be written as 



1 



N/2 

^ |m) <E> 

m=-N/2 



(4) 



Here |mi) <g> \rri2) is the eigenstate of J+ with eigenvalue 
mi + ui2 • The state vector \tjjme) is the only simultane- 



ous eigenvector of J+, J~ and J~ with null eigenvalue 

pera- 
and, 



ous eigenvec 

El El El- Obviously" the variances of these opera- 



tors vanish in the maximally entangled state \ipr, 
moreover, quite trivially, 

(Jz )me — ( Jy )me — i^x )trte — 7 



(^) 



where, here and henceforth (^4)^ = Tr[Agj], where gj is 
a generic density matrix, and g mt 



\l/jme)(lpme\ is the 
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density matrix of the maximally entangled state. Eq. (JjjJ 
provides therefore a necessary (but obviously not suffi- 
cient) condition for the realization of the maximally en- 
tangled state and will be exploited in the following to set 
up the feedback scheme. 

The QND measurement of the total population dif- 
ference between the two lower atomic states by photo- 
detection corresponds to the measurement of the observ- 
able J+ (the explicit form of the QND interaction Hamil- 
tonian is given by Eq. JHJl below), and can be realized, 
for example, by using light with two polarization compo- 
nents that interact differently with the two atomic states, 
generating in this way different phase shifts that produce 
a polarization rotation signal. Another method may con- 
sist in using modulated light, with one frequency compo- 
nent closer to resonance than the other, so that the inter- 
action with the atoms yields a phase difference between 
the two components. However, for ease of presentation, 
we follow the schematic model described in Ref. [TT1 | . 

The two atomic ensembles are placed in one arm of 
an interferometric setup, cf. Fig. The incoming field, 
whose annihilation operator is denoted by d+ , is a highly 
collimated single photon pulse, which is decomposed in 
two components by means of a 50% — 50% beam splitter: 
the reflected component ap follows the free path, while 
the transmitted component d/ goes through the atomic 
samples and interacts with them. 

Since only the state \g) is off-resonantly coupled with 
the excited level by the interacting field, the phase shift 
between the two field components is proportional to 
the population difference mi + mi between the atomic 
ground states and can be resolved by the intensi- 
ties measured in the two output ports of the interferom- 
eter by photo-detector D + , absorbing photons of mode 
d + = (ap + ai)/y/2, and photo-detector D_, absorbing 
photons of mode d_ = (ap — d/)/\/2. The sequence 
of measurements of the field phase shift corresponds to 
QND measurements of the observable J+ and yields a 
nondestructive evolution of the global state of the two 
atomic ensembles [Tlj . 

In order to guarantee a high spectral and directional 
resolution of the entering photon, we can place in each 
arm of the interferometer two symmetric ring cavities, in 
one of which the two atomic samples are enclosed (see 
Fig. 2). The ring cavity in the "empty" arm of the inter- 
ferometer is needed to achieve optimal mode matching 
at the output beam splitter of the interferometer. This, 
together with the condition of a detuning A much larger 
than the excited level decay rate 7, greatly reduces the 
effect of spontaneous emission which can then be omitted 
from our analysis. Moreover, by considering a traveling- 
wave probe in this far-off resonance case, any information 
about the relative positions of the atoms of the ensembles 
due to recoil is erased. Hence, the effective Hamiltonian 
of the atom-photon interaction can be written as |Tlj 

H ef f = ^ K (N-J+)a\a I , (6) 



where g is the coupling strength between the single atom 
and the radiation field. 

The modification of the state of the two samples in- 
duced by the photon-atom interaction and the subse- 
quent QND measurement can be formalized by the ac- 
tion of an appropriate POVM [3. Following Ref. [n[, 
we assume that photodetectors D+ and D- are able 
to reveal one photon at a time. Therefore, the state 
of the two atomic ensembles after each photo-detection 
is determined by the action of the (Kraus) operators 
M± = (I±e- lx( - N ^ J ^)/2, where I is the identity opera- 
tor and \ = g 2 r/ A is the phase shift, (t is the duration 
of the photon-atomic ensembles interaction process). 

The dynamics induced by the POVM, which preserves 
the rotational symmetry of the atomic system, and hence 
the value of the total angular momentum, cannot produce 
an evolution towards a maximally entangled state, which 
is a linear combination of eigenstates of different total 
spin J = Ji + J2 . The rotational symmetry can be broken 
by rotating, with opposite angles and — Q, each atomic 
sample around the a;— axes in the spin space. Actually, 
such a rotation is currently used in experiments for purely 
practical reasons, for example to clean up the relevant 
signal from technical noise at high frequency ||. The 
rotation is performed after each photo-detection, and is 
realized by the operator Ur = exp[— ifl J x \] exp^SlJ^] = 
exp[— i£lJ~]. In this way what is effectively measured at 
the n-th photo-detection is the rotated operator 

J+(nO) = e inQj * J+e- inQj * 
= (Jzi + J Z 2)cos(nn) 

+(J y l - Jy2)sm(nfl) , (7) 

as it appears in the non-rotated frame. 

Photo-detection losses are accounted for by consider- 
ing a finite efficiency of the measurement process, i.e. 
assuming that only a fraction 77 < 1 of the probe photons 
is actually detected. In this non-ideal situation the evo- 
lution of the density matrix is timed by the rate at which 
the single photon enters the interferometric set-up, and 
is conditioned by the possibility of photo-detection. 

If Q n is the density matrix of the total atomic spin sys- 
tem after n photons were sent on it, the (non-normalized) 
state g' n+ i at the successive step is given by 

Q' n+ i = < M±QnMl^ (8) 

if the photon is detected by detector D + (M+), or D_ 
(M_), with probability P± = r]Tr[M±g n Ml.} respec- 
tively; otherwise, if, with probability 1 — 77, no photo- 
detection occurs, it can be written as 

&+i = e~ iaJ ^ (M+QnMl + M-Q n Ml) e inJ * . (9) 

III. THE FEEDBACK SCHEME 

Our purpose is to efficiently realize the maximally en- 
tangled state g me in a controlled fashion. As shown in 
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Ref. ^lj, the measurement scheme described in section 
m yields only a certain probability that the global state 
of the two atomic ensembles will be gradually projected 
onto g me . Moreover, we recall that, to this aim, Eq. (JSJ 
is only a necessary condition and, although in the initial 
state go we have (J£)o = (Jy}o = (Jx)o = 0, as the 
measurements process goes on, the state of the two sam- 
ples does not necessarily satisfy this condition any longer. 
In fact, since we are measuring J+, the back-action of the 
POVM has the effect of decreasing the uncertainty of J+ 
and, at the same time, randomly shifting from its 

initial null value [HI El- As shown in Ref. Q, this 
shift can be interpreted stochastic rotation of the 
mean collective atomic spin around the y-axis induced 
by the measurement process. The optimal situation we 



would like to require is (J 2 f ) n = (J v 



at 



each step of the measurement process. We choose this 
requirement because it forces the state of the two atomic 
ensembles to satisfy condition JSJ at each stage of its con- 
ditional evolution and, due to the fact that \ipme) is the 
only simultaneous eigenvector of , J~ and J^T with 
null eigenvalue, this procedure likely increases the prob- 
ability to drive the system into the maximally entangled 
state. In fact, the validity of this choice is verified a 
posteriori by observing that the variances of the three 
collective operators J+, Jy and J~ arc monotonically 
decreasing functions during the time evolution. There- 
fore, to implement the optimal condition, Eq. ©, we im- 
pose a quantum feedback that properly counter-rotates 



the atomic spin operators. From the expression J7J) of 
the rotated operator at n-th step and photo-detection, 
it follows that the optimal condition on and J~ is 

obtained by imposing ( J+(nf2)}„ = 0. The simplest way 
to achieve this situation is to act with a unitary feedback 
operator Uf = exp [iX Jy] 0|. This operator realizes a 
rotation of an angle A around the instantaneous y— axes 
in the non-rotated frame and acts simultaneously on both 
atomic ensembles after the photon has been detected. 
This rotation can be realized, for example, by applying 
the combination of a static and an amplitude controlled 
rf magnetic field, which couples the two ground states 
\g) and |/) and drives the y-components of the collective 
spin operators generating the rotation Uf = exp [iX J+] 
[l4|. However, in our scheme, the amplitude of the driv- 
ing field is controlled by the feedback parameter A, that 
depends upon the measured signal. The state of the sys- 
tem after the photo-detection, the rotation around the 
"fixed" x-axes, and the action of the feedback, is given 
by the controlling updating equation 



Qn+l 



e i\J+ e -iUJ- 



M± Qn M\e 



(10) 



Obviously, if no detection occurs any feedback signal is 
implemented, and the updated state is given by Eq. 10. 



0, we get the 



Imposing the condition Tr 

required form of the feedback parameter A (A±), ex- 
pressed by the relation: 



{Jt)n + (e-^(«)j+ e *x4 + (n))„ ± (j+ e -ix[N-J+(n)] +e ix[N-J+(n)]j+} n 



The phase shift x ' 1S typically very small and when 
\N <C 1 (which is easily satisfied for N ~ 100), we can 
expand the numerator and the denominator of tan A± 
to second order in \. Since in the non-rotated frame 
the action of the feedback imposes (J+(nf2))„ = at 
each step of the photo-detection, Eq. Q implies that 
{Jt) n = (Jy)n = 0, and, as long as our scheme holds, 



is naturally satisfied. Retaining the leading 



terms in xN : we obtain the explicit expressions of the 
feedback parameter for the two possible photo-detection 
outcomes: 



tan Aj 



-X 2 iV(2((J+) 2 )»cosr!+([J-,J+] + ) ra sin») 
4(Jx + )„ 

2((J+) 2 )„cosr!+ ([J-,J+] + ) nS infi 
tanA_ = . - , . • (12) 



N(J+ 



with [A,B} + = AB + BA. 



IV. NUMERICAL SIMULATION AND THE 
ADIABATIC FEEDBACK 



We can now test by numerical simulations the measure- 
ment and feedback schemes. We simulate the conditional 
evolution of the global state of the two atomic samples, 
due to the acquisition of information by successive single 
photodctections, according to the updating equations for 
Qn+i, where the angle of the feedback rotation is deter- 
mined step by step either by the first or by the second of 
Eqs. Ifllty. 

For each case we have simulated 50 evolutions, with 
5 x 10 4 -j- 10 5 incoming photons. For every simulation 
we compute the overlap with the maximally entangled 
state Tr[g n g me ] and evaluate the degree of entanglement 
between the two atomic ensembles. Since in the case of 
perfect detection (rj = 1) we always deal with pure states, 
the entanglement is quantified by the von Neumann en- 
tropy E = —Tr[gj log 2 Qj] associated with the reduced 
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FIG. 3: (Color online). Numerical simulations of the evolu- 
tion in the presence of feedback of the overlap (Tr[g„g roe j) 
between the state of the two atomic samples (TV = 10) and 
the maximal entangled state in a case in which the feedback 
scheme is successful (S, dotted line) in generating g me and 
in a case in which the scheme is broken (73, solid line) by 
"anomalous" feedback angles (see text); inset: Evolution of 
the entanglement E in the S case (dotted line) and in the B 
case (solid line). The values of the parameters used in the 
simulations are: r\ = 1, x = 0-03 (achievable with highly col- 
limated beams), Q — n/10 and the total number of incoming 
photons is n p h = 5 x 10 4 . 

density matrices of the two subsystems of N atoms, Qj. 
If we restrict ourselves to states which are symmetric un- 
der permutations inside each subcnscmblc, E takes values 
between zero for a product state, and log 2 (iV + 1) for a 
maximally entangled state of the two samples. 

If 77 7^ 1 (mixed states) the entanglement can be ef- 
ficiently quantified by the relative violation of the lo- 
cal uncertainty relations (LUR) for the relevant ob- 
servables Jt, Jy and J~ , as shown by Hofmann and 
Takcuchi [ljj for entanglement generation in (N + 1)- 
level systems. This quantity is defined as Clur = 
1 — (5 J+ + 5 J~ + 8J~ ) /N, where 5 A denotes the variance 

of A and Clur. = 1 for maximally entangled states. 

We quantify the efficiency of the scheme by determin- 
ing the fraction of simulations T for which the value of 
the overlap between the final state of the two atomic en- 
sembles after the photo-detection sequence and the max- 
imally entangled state is larger than 0.99. The results of 
the simulations show a substantial efficiency of the feed- 
back scheme (FS) compared with those obtained by pure 
probabilistic schemes (PS) without feedback. We now 
discuss the different cases: 

• a) No feedback: in this (PS) case the successful rate 
is very low: in fact, we have T = 6% both in the 
ideal case (rj = 1) and in the more realistic case 
(77 = 0.9). 

• b) Simple feedback: the action of the simple fecd- 



FIG. 4: (Color online), (a): Numerical simulations of the 
evolution of ( (solid line) and of 8 j£ + SJy (dotted line) 
in the S case; inset (a): evolution of the feedback angle A 
during the S simulation; (b): Evolution of (J^) (solid line) 
and of <5Jjf + 8 Jy (dotted line) in the B case. Notice that 
after the point in which (J^) < 8jf + 8Jy, the variances 
explodes. This corresponds to exploding values of A around 
that point (inset (b)). Parameter values are the same as in 
Fig. 3. 
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FIG. 5: (Color online). Evolution of Clur in a S simulation 
(dotted line) and in a B simulation (solid line). Positive values 
of Clur provide a quantitative estimate of the amount of 
entanglement. Parameter values are the same as in Fig. 3 
except that now 77 = 0.9. 



back relevantly increases the successful rate with 
respect to the pure PS. In fact, in the ideal case of 
perfect detection we have T sa 54%: the action of 
the feedback on the mean value and the variances 
of the relevant spin operators j+, J~ and J~ is 
extremely effective, driving them to the ideal null 
values as the state of the two atomic clouds (g n ) ap- 
proaches to the maximally entangled state (g me ), 
cf. Fig. ^a) and^Jb). In the more realistic case 
(rj = 0.9) we obtain T « 24%. The worse result 
in this latter case is obviously due to the the fact 
that the instability of the scheme is even stressed 
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FIG. 6: (Color online), (a): Evolution of the overlap in the S 
case (dotted line) and in the B case (solid line); (b): Evolution 
of the purity Tr[g^] in the S case (dotted line) and in the B 
case (solid line). Notice that a successful feedback scheme 
has also a purification action. Obviously when the scheme is 
broken the purification effect is lost. Parameter values are the 
same as in Fig. 5. 



by the mixedness of the state of the two atomic 
ensembles, making the action of the feedback less 
"precise" than the case with pure state. Moreover, 
the evolution of the entanglement is slower, requir- 
ing many more photons to be sent on the samples 
to get the same values of the overlap of the case 
with rj = 1, cf. Figs. |5] and It is remarkable 
that, since the feedback signal depends only on the 
statistical moments of the spin observables, its ef- 
ficiency does not depend on the exact values taken 
by the feedback angles. In fact, extended numeri- 
cal simulations show that the feedback mechanism 
is not appreciably affected by perturbing the values 
of A by fluctuations of the order of A± x 10 _1 H-10 -2 . 
This is an important aspect of our scheme, which is 
therefore quite robust against imperfections in the 
actuation of the feedback operation. 

c) Modified ( "adiabatic" ) feedback: The only par- 
tial success of the feedback scheme is caused by the 
possibility that Eqs. (|12|) provide wrong values of 
the feedback angle when g n w Q me - In fact, in this 
situation, the feedback angle A should uniformly go 
to zero, otherwise, since g me is not invariant under 
the feedback rotation, the state of the two sam- 
ples is transformed in a state practically orthog- 
onal to the maximally entangled one, cf. Fig. 
This case could happen when ( J+) becomes smaller 
than the numerators in Eqs. i|12|) . essentially pro- 
portional to (5J+ + S J~, cf. Fig. 0](b). Actually, 
from the uncertainty relations for spin systems is 



to show that 5J+SJ,, 







(Ji) 







and, due to the random jump caused by 
the back action of the photo-detection, in this sit- 
uation it is possible to have (J+) < SJ+ + SJ~, 
getting in this way "exploding" values for A, cf. 
inset Fig. 0] (b). However, it is possible to avoid 



these "anomalous" feedback angles by modifying 
the feedback algorithm in order to force A± to go 
uniformly to zero as the maximally entangled state 
is approached. In fact, inspired by the mean values 
taken by A± in simulations in which the maximal 
entangled state is reached (inset Fig^Ja)), we can 
impose on the feedback angles of Eqs. (|12|1 a de- 
creasing exponential cut, a sort of adiabatic switch- 
ing off, expressed as a function of the number n of 
incoming photons: more precisely, the exponential 
cut is X cu t = ircxp(— x), with x = n x 10~ 4 , and 
we impose that if |A±| > X cu t, then A = e(A±)A cll t, 
where e is the signum function; otherwise A = X±. 
However, in order to not reduce the effectiveness 
of the feedback scheme, we switch on the correc- 
tion procedure only after a high degree of entangle- 
ment is obtained, so that the state of the atomic en- 
sembles is "close" to the maximally entangled one. 
The point in which the correction procedure should 
begin to act has to be chosen in such a way that 
the probability that anomalous values of A have al- 
ready occurred is very low: by a statistical analysis 
of the simple feedback scheme, this point coincides 
roughly with n»2x 10 4 . In this way we can softly, 
even though more slowly, drive the state of the sys- 
tem toward the maximally entangled state. With 
this procedure, the results of the numerical simu- 
lations are remarkable. In the case with r\ = 1, we 
have T « 92%, and the successful rate remain re- 
markably high also for the realistic choice r\ = 0.9: 
T w 78%. The reliability of the scheme can also 
be appreciated by considering stronger constrained 
data: in fact, if we consider overlaps greater than 
0.999, we have 90% for ^ = 1 and 64% for rj = 0.9, 
while overlaps greater than 0.9999 (practically, a 
full realization of the maximally entangled state) 
leads to 80% for -q = 1 and 50% for -q = 0.9. 



V. CONCLUSIONS 

In conclusion, we have introduced an experimen- 
tally feasible, quasi-deterministic stroboscopic feedback 
scheme able to generate the maximally entangled state 
of a system composed by two mesoscopic atomic ensem- 
bles, each containing a number of atoms N ~ 10 — 10 3 . 

The feedback scheme is based on QND measurements 
of the total atomic population difference between the in- 
ternal states of the two atomic samples, which are probed 
by a photon beam. Though the state of the atomic sys- 
tem is conditionally entangled by the photodetections, 
the conditional evolution docs not necessarily lead to the 
the maximally entangled state. To this aim, we designed 
a scheme in which the signal obtained by the measure- 
ments, and the corresponding information on the state 
of the ensembles, is exploited to set up a feedback signal, 
formalized by an unitary operator acting on the state of 
the atomic samples, that drives the system into the max- 
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imal entangled state. The feedback strategy adopted is a 
modification of that introduced in Ref. [lj| by Thomsen, 
Mancini and Wiseman. 

The procedure is analyzed by numerical simulations of 
the conditional evolution of the state of the two atomic 
samples. In particular, for every simulation, we have 
monitored the overlap between the state of the atomic 
system and the maximally entangled state, and we eval- 
uated the degree of entanglement between the two atomic 
ensembles. The efficiency of the protocol is quantified by 
the fraction of simulations in which the final value of the 
overlap is larger than 0.99. 

The results of the simulations show that the imple- 
mentation of the feedback scheme yields a considerable 
enhancement of the efficiency of generating the maximal 
entangled state of the two atomic ensembles, compared 
with the case of pure probabilistic schemes without feed- 
back. Moreover, we have shown that this efficiency can be 
further improved by modifying the procedure by means 
of an adiabatic switching off of the feedback signal. 



We have shown that the scheme is able to tolerate a 
number of experimental imperfections, as, for instance, 
a limited detection efficiency. In the numerical simu- 
lations of the feedback protocol, we have considered a 
photocounting efficiency r\ = 0.9: although this value 
is optimistic (typical laboratory photocounters are char- 
acterized by r\ ~ 0.5), the impressive, continuing pro- 
gresses in this field give hope that very high efficien- 
cies will be available in the near future (see for instance 
Ref. [2(j, where r\ ~ 0.88 has been achieved at a par- 
ticular wavelength). Nevertheless, at variance with the 
experimentally achieved 0| probabilistic generation of 
weakly entangled states between macroscopic atomic en- 
sembles (with typical numbers of atoms N ~ 10 13 ), our 
quasi-deterministic scheme applied to mesoscopic ensem- 
bles allows to access a presently unexplored area of entan- 
glement production in atomic systems, even with present- 
day technologies against photon losses and imperfect de- 
tections. 
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